On knotted streamtubes in incompressible hydro dynamical flow and a restricted 

conserved quantity 
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For certain families of fluid flow, a new conserved quantity - stream-helicity - has been established. 
Using examples of linked and knotted streamtubes, it has been shown that stream-helicity does, in 
certain cases, entertain itself with a very precise topological meaning viz., measure of the degree 
of knottedness or linkage of streamtubes. As a consequence, stream-helicity emerges as a robust 
topological invariant. 
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I. INTRODUCTION 



II. STREAM-HELICITY 
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Lord Kelvin (who alongwith Hclmholtz pioneered in 
the subject of vortex motion) had recognised in the late 
19th century that in an inviscid and barotropic fluid be- 
ing acted upon by irrotational body forces, any linkage or 
any knottedness in the vorticity field at any earlier time 
should remain conserved at all later times. After almost 
hundred years, MoreaufH and later Moffatt Q established 
an invariant known as helicity which is of topological 
character and encompasses Kelvin's insight. Stark anal- 
ogy between vorticity (c3) in ordinary fluid dynamics and 
magnetic field (B) in magneto-hydrodynamics (MHD) 
prompted Moffatt to give similar topological interpreta- 
tions to magnetic helicity and cross-helicity (which, by 
the way, measures the degree of 'mutual' knottedness of 
two the fields: w and B). Hence, researchers were able 
to effectively connect the two very rich fields viz., topol- 
ogy and fluid dynamics and excited a lot of interest in 
this direction. But what lord Kelvin had missed was the 
possible existence of knotted streamtubes in the steady 
Euler flows, a fact very logically speculated by Moffatt Q. 
Not much has been done on that. Here, in this paper, in- 
spired by the works of Moffatt, we shall introduce in sec- 
tion (II) a quantity, which we shall call 'stream-helicity' 
(S) in inviscid and incompressible fluid being forced by 
irrotational body forces. It will be shown that stream- 
helicity is a conserved quantity under certain restrictions 
which are not, of course, very rare in practice. In section 
(III), we shall note how this conserved quantity can have 
a very sound topological meaning for at least some kinds 
of flows and hence, how stream-helicity can be raised to 
the status of a topological invariant for linked and knot- 
ted streamtubes. 



Let us start with the Euler equation (equation (fTj)) 
for three-dimensional, inviscid, incompressible fluid being 
acted upon by irrotational body forces. P used in the 
equation includes the effect of such forces also. Since the 
fluid is incompressible, i.e., the density is constant, we 
are taking the density to be unity for convenience. 
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u + (u.V)u 
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Incompressibility yields for the velocity field u: 
V.u = 



(1) 



(2) 



which helps in defining the vector potential £ for the ve- 
locity field as follows: 



V x £ 



(3) 



Obviously, £ is not unique, for, a term VA, A being 
a scalar field can always be added to it keeping u un- 
changed. We shall come back to this issue in the right 
place. For now, let us put relation ([3]) in the equation 
© to get: 



d 



dt 

But we have: 
[V x (tf.V)l], 



(V x + («.V)(V x i) 



-VP 



(4) 



Using relation 
d 



= €ijkdj(uidi£ k ) 

= e ljk (djUi)(di^ k ) + tijkUidjdi^k 

= ey fc (S,-u,)($&) + [(tT.V)(V x f)]^) 

in the equation (j4|) we get: 



dt 



(Vx£) + Vx{(u.VK} 



= eijk(djUi)(di£k) ~ diP 
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=> -7^£ + (tt.V)£ = curl 77 
where if is defined as: 

m = eijk(djUi){diCk) - diP 



(6) 



(7) 
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Now, let us define ' stream- helicity' (S) as: 
S= [ luSx 



(8) 



where V is a volume occupied by the fluid. At this point 
let us ponder over the aforementioned non-uniqueness of 
the vector potential^. For smooth discussion's sake, we 
assume for the time being that the volume is simply con- 
nected. Suppose £j — > £j + diX, then from the definition 
([8|) of stream-helicity we can find the change SS in S to 
be: 



SS = / VX.ud s x = 



Xu.hd x 



(9) 



BY 



where n is the unit vector perpendicular to the infinitesi- 
mal surface element d 2 x and we have used the relation ([2]) 
and Gauss divergence theorem. The relation (J9j) amounts 
to saying that the stream-helicity will be gauge invariant 
in case the surface dV bounding V, is the surface made 
up of streamlines i.e., u.n = on dV . This condition for 
gauge invariance is rather strong because if u.n ^ on 
dV then one cannot seek refuge in Coulomb gauge for it 
is too loosely defined inside V with no information about 
the outside field whatsoever. More starkly, it means that 
different solenoidal vector potentials inside V can corre- 
spond to Coulomb potentials of fields which have differ- 
ent structures outside V. Now if we relax the condition 
that V is simply connected, then the line integrals of £ 
about the 'holes' in the possibly multiply connected re- 
gion have to be specified in order to have gauge-invariant 
stream-helicity within dV on which u.n = 0. The upshot 
is that only in the gauge u.n = (which in the simple 
language means that the boundary is impenetrable), S 
can be conserved. 

We now wish to demonstrate that under certain restric- 
tions this quantity is in fact conserved. So, we take total 
derivative of S w.r.t. time to get: 



dS 
~dt 
dS 
~dt 



^.u)d 3 x 
C{-VP)d 3 x 



u.(cmr 1 ff)d 3 x (10) 



where, D/Dt is the material derivative w.r.t. time and 
it basically is a shorthand for d/dt + u.V. Again, simple 
vector algebra suggests: 

(V x |).(curl -1 r?) = V.(|* x curl -1 17) + 7 7-£ ( n ) 

With relation ([3]) in mind, inserting relation (|TTj) in the 
equation (fTU|) . we have the following: 



2 J £.VPd 3 x + J V.(£ x curr 1 ff)d 3 x 
€idi(e ijk € k djUi)d a x (12) 



dS_ 

dt 



where, equation © has been used. The first two terms 
of the equation (fTSj) can be changed to integration over 



the surface which bounds the volume V in question (the 
surface will obviously extend to infinity if the fluid is 
unbounded) using Gauss divergence theorem and so if £ 
decays fast enough to go to zero on the bounding surface 
then these two term vanish; there may be other reasons 
for the terms to vanish as will be seen in the next section. 
Now, let us consider the third term. If this term vanishes 
then only one may set 



dS 
~dt 



= 



(13) 



and say that stream-helicity is a conserved quantity. 
Though it seems to be very restrictive, but one can see 
that in the following commonly occurring cases the inte- 
grand of this term trivially vanishes. 

1. The vector potential is one dimensional. 

2. The vector potential has no dependence on the di- 
rection along the velocity field. (Other conditions 
given below are basically this condition's corollary.) 

(a) £ is two dimensional but has dependence only 
on the third direction. 

(b) £ is two dimensional with spatial variations 
only on the plane containing it. 

(c) £ is three dimensional but depends only on 
any one of the three independent directions. 

One can see that such flows are very commonly found 
in any elementary text-books on fluid mechanics. For 
example, in accordance with case (1), for the one dimen- 
sional vector potential: £ = kfl(x 2 + y 2 )/^, the corre- 
sponding flow is u — —ifly/2 + jflx/2 which basically 
is the velocity field for three dimensional fluid counter- 
clock-wisely rotating about z-axis. As another instance, 
this time to go with the case (2) (2(b), to be precise), is 
that of a uniform flow along x-direction: u — Ui which is 
generated by the vector potential £ = —jUz/2 + kUy/2 
that evidently is two dimensional with spatial variations 
only on the y-z plane containing it. A rather non-trival 
case (as an example of the case 2(a)) is for the flow: 
u = iU sinaz + jU cos az (a is a constant) for which the 
vector potential is £ = i(U /a) sinaz+j(U/a) cosaz; read- 
ers must have noticed that this is just a variant of the 
more general flow viz., ABC flow (see e.g., 
So, for the families of fluid flow for which the vector po- 
tential falls into the above set and if eventually the equa- 
tion (TIU)) holds, stream-helicity is a conserved quantity. 
Also, for the fluid flows for which doesn't fall in the above 
set but the integration goes to zero some reason or the 
other (which has not been investigated), S will remain 
conserved. 



3 




by the linked structure only, then 



FIG. 1: Linked closed streamtubes. The tube and hence the 
streamlines inside are not twisted i.e., the fluid inside the 
tube doesn't swirl. The directions of arrows are showing the 
direction of the streamlines filling the tubes. 



III. TOPOLOGICAL MEANING OF 
STREAM-HELICITY 



Now, we ask the question if it is possible to give stream- 
helicity a topological meaning and more importantly, can 
that topological meaning turn out to be a topological in- 
variant. We shall see that the answer is in affirmative. 
To get both of the expectations met, one (other uninves- 
tigated possibilities may also be there) of the ways seems 
to be the following: 

Consider two circular thin streamtubes which are singly- 
linked and for the two tubes the 'strengths' are Vc x and 
Vc 2 respectively, where C\ and C2 denote axis-circles of 
the corresponding tubes. By 'strengths' we mean that 
Vc x = u.ddi and Vc 2 = u.da.2- (See FIG-1). Again, we 
assume that the velocity field, we shall be dealing with, 
is generated by a vector potential £ which is a Beltrami 
field i.e., 



V x £ = a£ 



(14) 



where, a is a numerical constant. Moreover, suppose that 
of the conditions gathered in the previous section for £, 
at least one is applicable, say the second one that: 



(tf.V)£ = 



(15) 



whether this is possible or not may be a valid question. 
One may derive 'some' relief from the fact that if £ is anal- 
ogous to ABC flow (Gromeka(1881); Beltrami(1889)), 
then it does satisfy such condition though unfortunately 
it may not sustain a linked structure of streamtubes. 
Then the streamtubes will be made up of streamlines 
which are coincident with the 'flux-lines' of £ field. 
If we define the volume over which the integration is de- 
fined for the stream-helicity to be the volume occupied 



S 



S 



£.ud 3 x 



1st Strcamtubc 

£.ud 3 x 



V Cl / £.dli 



2nd Strcamtubc 

Vc 2 



Z.dl 2 



= 
= 



Vc t / / u.da 

J JDCx 

V Cl Vc 2 +V Cl V C2 
2V Cl V C2 



G 2 

Vc 2 



'.da 



DC 2 



(16) 



where in the preceding steps we have used ud 3 x — > 
Vc ± dli, Vc 2 dli on C\ and C2 respectively; DC\ and DC2 
denote the area spanned by C\ and C2 respectively. Ob- 
viously, if the linking number is n and not one as in this 
case, one would easily generalise the result to: 



S = 2nV Cl V C2 



(17) 



which, being dependent on the mutual linking of stream- 
tubes, is a topological quantity. One may write from 
equation (fT2")l using Gauss divergence theorem in the fol- 
lowing form: 



dS_ 
~dt 



-2 I (Cn)Pd 2 x+ I (curl _1 ^.(n x ^)d 2 x 
£idi(ei jk £ k djUi)d 3 x (18) 



where fi is the unit vector perpendicular to the surface 
at each point on the surface of the linked structure. The 
first term and the third terms of the equation (118|) are 
zero in this case by construction of the linked structure; 
so is the second term but it needs a bit of manipulation 
as explained below. 

First of all, we use equation to rewrite the integrand 
of the second term of the relation (TT51) as 



at 



.(n x I) = e ijk nj£ k 



dt 



(19) 



Now, if we consider the Frenet-Serret coordinate system: 
(T,N,B), then in the case we are considering £/|£| = T 
and h = N; obviously on the surface of the specific tube 
we are considering, at each point, the triad so that there 
£2 = £3 = and ni = n 3 — and hence due to the 
antisymmetry of €ij k we have : 

(curr 1 ^).^ x I) = (20) 
So, obviously we land up on the following relation: 

dS , . 

-dt=° < 21) 

Therefore, for incompressible, ideal and conservatively 
forced fluid flow, in certain configurations, we can have a 
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topological invariant - stream-helicity - for linked struc- 
tures of streamtubes. 

So far so good. So, stream-helicity does seem to make 
physical sense for linked two (or more) streamtubes. But 
what if a single streamtube is knotted? A single knotted 
streamtube must have an unavoidable twist of velocity 
field (which we hope, in this case also, may be derived 
from a Beltrami velocity vector potential and is of similar 
kind as has been dealt with earlier in this paper) within 
the tube. How to deal with such a scenario has been 
discussed for knotted vortex filaments by Moffatt [(|. We 
know when an arbitrary tame knot is viewed in a stan- 
dard plane of projection with finite number of crossings, 
each of which is either positive or negative, it can be 
changed to a unknot (and ergo, subsequently continu- 
ously deformed to a circle) by switching the crossings 
for a finite number of times. (To remind the readers, a 
crossing is defined as positive or negative according as the 
overpass must be rotated counter-clockwise or clockwise 
to bring it into coincidence with the underpass.) One 
may note that the resulting circle may be converted back 
to the original knot simple by performing the operations 
in the reverse order. With this in mind, let us consider a 
tubular region with the circle as axis. The cross-section 
of the tube is small and over that the velocity of field, 
which we suppose is filling the tube with strength V , is 
uniform; each streamline is, of course, a concentric circle 
to the circle serving as the axis. Now, let us transversely 
cut the tube somewhere and reconnect it back after giv- 
ing it a twist through an angle 2tvN (where, N is an 
integer). This way we are introducing a stream-helicity 
of magnitude NV 2 . Then by introducing proper switch- 
ing loops with similar strength, this construction may be 
changed to a knot with stream-helicity: 

S= [N + 2(n+ -n_)]V 2 (22) 



where, n + and n_ are respectively the number of positive 
and negative switches needed to create the knot whose 
stream-helicity we are interested in. One may prove that 
N + 2(rt + — rt_) is actually the linking number of any 
pair of streamlines in the knotted streamtube. It also 
is the self-linking number of the 'framed' knot which is 
'framed' using Frenet-Serret coordinate system [7]]. Point 
to be noted is that for the kind of velocity field we are 
discussing S will remain conserved and hence emerges as 
a topological invariant, for, evidently S depends on the 
topology of the knotted streamtube. 



IV. CONCLUSION 

To summarise, a new conserved quantity - stream- 
helicity - has been found (albeit, in a restricted sense) in 
fluid dynamics. By seeking a topological interpretation 
for it in the certain configurations of linked and knotted 
streamtubes, the bridge between topology and fluid 
dynamics has been made even stronger. In addition, 
as a byproduct, the seemingly non-physical quantity 
- velocity vector potential - has given itself a sort of 
physical meaning by getting involved in measuring the 
degree of knottedness of streamtubes. 
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